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Abstract 

R parity violating trilinear couplings of the minimal supersymmetric 
standard model (MSSM) are constrained from the limit on the electron neu¬ 
trino mass. Strong limits on these couplings follow from the earlier neglected 
contribution due to sneutrino vacuum expectation values. The limits on most 
of the couplings derived here are stronger than the existing ones for a 
wide range in parameters of MSSM. These limits strongly constrain the in¬ 
terpretation of recent HERA results in e'^p scattering in terms of production 
of squarks through R violating couplings. In particular, the interpretation in 
terms of production off strange quark as suggested recently is not viable 
for wide ranges in MSSM parameters. 


The baryon and the lepton nnmber symmetries are known to be violated when the standard 
electroweak model is extended to inclnde snpersymmetry. These violations arise throngh 
the following terms in the snperpotential of the minimal snpersymmetric standard model 
(MSSM): 


IT^ = -Ah, L'Q'D- - A 


If 
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u^D'; 
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( 1 ) 


where Qb, and L'^ denote the snperhelds corresponding to the weak eigenstates of qnarks 
and leptons respectively. The presence of the lepton nnmber violating interactions leads to 
interesting deviations from the standard model snch as —>■ vr+i/z/, nentrino masses, 

violation of weak nniversality etc. The experimental limits on these have been used to 
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constrain [Q] the relevant parameters in eq.(l). A particnlarly stringent constraint follows 
from the limit on the electron nentrino mass. The existing analysis |^,^] of these constraints 
has neglected an important and dominant contribntion to the nentrino masses which wonld 
occnr when any of the lepton nnmber violating term in eq.(l) is present. Incorporation of 
this additional contribntion can change 0] the already obtained limits significantly. The 
aim of this paper is to systematically derive these constraints and discnss its implications. 

A carefnl analysis of the constraints on X' conplings is of particnlar signihcance in view 
of the recent anomaly observed by the HERA experiments HI P] and ZEUS [||. They have 
reported excess events in deep inelastic e~^p scattering at > 1.5 x 10^ GeV^. This excess 
cannot be reconciled with the expectations based on the QCD improved parton model and 
point to the existence of physics beyond standard model. The reported excess conld be at- 
tribnted to the presence of non-standard contact interactions 0, to the resonant prodnction 
of a leptoquark state ^ or of a sqnark throngh its R parity violating conplings A' PJT0[|. 

The possible resonance seen at HERA can be identihed either with a or a sqnark of 
the MSSM. Three plansible prodnction mechanisms have been proposed and discnssed: 
e^dR —>• Cl, e.'RdR tL, e.'RSR ii- All these mechanisms (particnlarly the last one) reqnire 
signihcantly large valnes for which are constrained among other things by the electron 
nentrino mass. 

The nentrino masses arise throngh two different sonrces [|TI| in the snpersymmetric stan¬ 
dard model. The lepton nnmber violating trilinear terms in eq.(l) constitnte the hrst sonrce. 
They generate majorana masses at the one-loop level ||TI|JT^ . The lepton nnmber violating 
terms in snperpotential are invariably accompanied by similar terms in the scalar potential 
providing another sonrce of nentrino masses. Specifically, the following terms are allowed 
by SU{3) X SU{2) x f/(l) invariance in V: 


Usott — —Bj}. Uild^ + rnt^TT, + c.c -!-••■ 


( 2 ) 


These terms are present at the GUT scale if original snperpotential contains bilinear R 
violating terms. Bnt terms in eq. (2) wonld be generated P,|I^ at the weak scale even when 
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the violation of lepton number is only through the trilinear terms as in eq.(l). The major 
consequence of these terms is the generation of sneutrino VEV which causes the mixing of 
neutrinos with gauginos leading to the neutrino masses P,p!3|-[T6[ . 


The bounds on A'lp quoted in the literature are derived using the loop induced 
contribution to the neutrino masses. Recent analysis of neutrino masses P,P!B|-[1T[] has re¬ 
vealed that the neutrino gaugino mixing tends to generate a larger neutrino mass compared 
to the loop induced contribution. This contribution can thus lead to stronger bounds on 
the trilinear couplings. We derive these bounds. It turns out that for a range of the MSSM 
parameters the bounds on following from the neutrino mass can be stronger than the 


existing bounds based on other processes such as —>• vr+z/n |]^, neutrinoless double beta 


decay Z —^ e'''e [IT^ etc. As a result, a sizable part of the parameter space otherwise 


allowed for the interpretation of HERA events is curtailed by the neutrino mass bounds. 

We work with R violating version of the MSSM with a standard set of soft supersymmetry 
breaking terms specihed at a high scale near Mcut = 3 x 10^® GeV. Since our aim is to derive 
bounds on A' couplings, we concentrate only on that term in eq.(l) and set other couplings 
to zero. It is more convenient to express this equation in terms of the mass eigenstates 
Qi, D'^ and of quarks and charged leptons respectively . 


IVn = -A' - Ah, , 


( 3 ) 


where Ah, are related to Xijk of eq.(l) and 


\u _ -i/KM \f 

\jk — ^Ij \lk 


The lepton number violating part of the scalar potential at the weak scale is given by 


V^oft = - B, + mlH, i>*H, + c.c + -- - , 


( 4 ) 


i^Hi 


( 5 ) 


where we have explicitly written only those terms which are responsible for generating the 
electron neutrino = z/) mass and Ah = The terms in eq.(5) generate a non-zero 
VEV(z/) and give rise to the following neutrino mass |Tl| J l3|-jl6 i ^ : 
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{cg^ + 9 '^) {!') 


i^^Ug'jtree — 


4 (—c fiM + COS (3 sin (3 {c + tan^ 6w)) 

where c ~ bg'^/?>g'^, g"^ and g'^ are gauge couplings and tan/3 = {H 2 )/{H^) = V 2 IV 1 . The 
(h) follows from the minimization of the scalar potential which contains terms of eq.(5) in 
addition to the standard soft terms of MSSM. Explicitly, 


(^) 


B^V2 - 

ml + cos 2(3 


(7) 


ml is the soft SUSY breaking sneutrino mass. 


All parameters in eq.(6) are specihed at a low scale Qq . We have chosen Qq to be Mz 
2l|] . The parameters Bi,, m^j^ are assumed to be zero at Mgut ■ Their values at Qo follow 


from the following renormalisation group(RG) equations 
dB 


dt 

dt 
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( 8 ) 


(9) 


In the above equations, t = 2 In^McuT/Qo)] hj?, hf are Yukawa couplings corresponding 
to the down,up quarks and the charged leptons respectively ; = (h^’^’^)^/167r^ ; 

dj = tti/dvr . The sum over generation indices fc, I is implied above. The parameters and 
satisfy the following RG equations: 


dX 


kk 


dt 

dA 


kk 


dt 


Kk {^^^k ~ 2^k + ^“1 + 2^2 + , 

3 9 1 16 7 

^^IkYk^ - -^^kYk - - y«3M3 - 3d2M2 - -a,M, . 


( 10 ) 


We have kept only leading order terms in A/!;, in writing eqs.(8 -10). Rest of the param¬ 
eters appearing in eqs. (8 -10) satisfy the standard RG equations to this order in Note 
that the second term in eqs.(8-9) generate non-zero By and ml^j^ at Qo- The neutrino mass 
i'knyjtree following from eqs.(6 -10) thus involves the combination {X^.j^hp^. 


The trilinear interactions in eq. (3) lead to the following my^ at the one-loop level ||TI|,p!2 
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{m^Jioop ^ sin 0 fccos 0 fc In-^ , 

lOTT iVi^^ 

where we have implicitly assumed that only one X[jj^ is non-zero at a time. <pk and M|^ 
respectively denote the mixing among squarks dk, and their masses. The mixing cpk is 
proportional to and hence the {niujioop also scales as therefore provides 

a bound on which can be converted to a bound on X'^^^ using eq.(4). The resulting 
bounds become stronger with increase in tan/? due to the fact that involve = 

/V cos j3 . For the same reason, bounds display strong hierarchy, typically 

{^kk)max n n\ 

fA'"') "" 

K'^lllmax ill'k 

It also follows from eq.(4) that 


(Kkk) max 
{^Uk)max 


■frKM 

^Ik 

\rKM ■ 
^kk 


It is seen from last two equations that the A'^ 33 , (A 333 ) is constrained most (least) by niu^. 

The Kobayashi Maskawa mixing neglected in the earlier analysis plays quite an important 
role here which needs to be emphasized. The earlier analysis could not constrain A^^^ with 
j ^ k due to the fact that single such A^^i. leads to zero when Vkm is neglected. This is 
no longer true in the presence of Vkm allowing us to derive signihcant bounds on couplings 
^ijk (i 7^ ^)- 

We have numerically solved eqs.(8-10) and the standard equations of other parameters 
appearing in them. They lead to and at to = 2 In which are specihed in terms 
of the standard MSSM parameters which we choose as gaugino (gravitino) mass, M 2 (m), 
tan/? and universal trilinear strength A. B^ and determine {m,y)tree in terms of these 

parameters and A^^(to)- The 1-loop contribution also gets hxed in terms of these parameters 
since <pk and ^ appearing in eq.(ll) are determined using the standard 2 x 2 mixing 
matrix for dk, d% system. 

The bounds on different couplings are displayed in Figs. 1 and 2 . Apart from being 
dependent on SUSY parameters, these bounds are quite sensitive to the chosen sign of 
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the /i parameter. This is due to the fact that for one (namely negative) sign of p, two 
terms appearing in the sneutrino VEV, eq.(7) cancel while they add for positive sign. The 
suppression in sneutrino VEV occurring in the hrst case weakens the bound on A'. Such 
suppression has been previously noted in the case of bilinear R violation P,p!l|-[T 6 [] . Fig. 
1 displays bounds on couplings and A']^ 2 i obtained by demanding = {m^Jtree + 
{jnv^)ioop < 5eV. Curves for three representative values of tan/3 and universal gravitino 
mass m are shown. For comparison, we also display in the same hgure the existing bounds 
on these couplings. The most stringent bound on is derived from neutrinoless double 
beta decay and on A '121 from the process K~^ —>■ 7t~^u9. These are shown as function of M 2 
in the same hgure using MSSM expressions for the relevant squark masses. It is seen that 
the bounds derived here are comparable or better (in case of larger values of M 2 ) than the 
already existing ones. 

Fig 2a (/X > 0) and 2b (p < 0) show the bounds on A '^33 and A']^ 32 . The A )33 is constrained 
most and is required to be as low as 3 x 10“® even for tan /3 = 5. For comparison we also show 
the limit on A '^33 which would follow if sneutrino VEV is completely neglected as in previous 
works [0,§|. It is clear that inclusion of this VEV drastically alters the bound. We also 
display limit on A 332 which was poorly constrained in the previous analysis and which 
is relevant for the interpretation of the HERA events. The limits displayed in £g 2a are at 
least an order of magnitude stronger than the value needed to reconcile HERA anomalies. 
The resulting bounds on other couplings A^^ not explicitly displayed in above hgures can 
be read off from eq.(13). These constraints are listed in a table for two representative values 
m = M 2 = SOGeV and m = 50 GeV, M 2 = 75 GeV and tan/3 = 10. These respectively 
correspond in MSSM to ~ 163, 240 GeV. It is seen that the neutrino mass considerably 
improves on the existing constraints. 
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\f 

'^ijk 

Previous limits 

m = 50 GeV 

M 2 = 50 GeV 

m = 50 GeV 

M 2 = 75 GeV 

111 

0.002 

18 


0.0012 

0.0052 

121 

0.024 

17 


0.0056 

0.023 

131 

0.52 [gg 


0.015 

0.065 

112 

0.024 

17 


0.0003 

0.0012 

122 

0.024 

17 


6.3 10-5 

2.6 10 -^ 

132 

0.56 1 

3 


0.0015 

0.0065 

113 

0.024 0 

2.3 10-5 

9.8 10-5 

123 

0.024 

17 


4.6 10-5 

1.3 10-^ 

133 

0.002 [§ 


1.9 10-5 

7.9 10-5 


Table 1. Limits on single \Lf. following from the electron neutrino mass for A = 0, 
tari/S = 10 and /i > 0. The limits become stronger for larger tan f3 and weaker for 
negative /i. The existing limits mentioned in the table are for the relevant squark 
mass ~ 200 GeV and gluino mass M 3 = 200 GeV in case of 


We now briefly discuss the implication of the above bounds on the interpretation of HERA 
events deferring details to a latter work. The HI and ZEUS data allow interpretation as 
resonance production of squarks provided P,pT 


('^121) '^131) ^1/2 ’ 


(14) 


where B refers to the branching ratio for the squark decay to ge"*" and f is the t^ — tn 
mixing angle. Limits on couplings displayed in Figs 1 and 2 are indeed stronger than above 
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(even for i? ~ 1) for a large range of MSSM parameters. In particular, the interpretation 
in terms of the production off strange quark at HERA is severely constrained from 
This possibility was studied in detail in which included constraints on B from CDF, 
LEP constraints on 6p and restrictions on A '^32 coming from Z partial widths. The last one 
gives fairly weak constraint on A 332 compared to the one coming from This constraint 
therefore places strong restrictions on the available parameter space. These are displayed in 
Fig. 3 where we show the region in m — M 2 plane in which neutrino mass bound is satished 
and A 332 > 0.3. We have also displayed the restrictions coming from the experimental bound 
of 85 GeV on the chargino masses at LEP. The combined allowed region becomes smaller as 


tan (3 increases and eventually beyond tan /3 = 10 it gets pushed to m > 240 GeV This 
will be limited further when one includes other constraints such as Sp and GDF constraints 
on B relevant, when the squark mass is below 210 GeV. The corresponding restrictions are 
much stronger when p > 0. For example, we hnd in this case that neutrino mass bound is 
inconsistent with A 333 > 0.3 for m, M 2 < 600 GeV and tan/? > 1 if A = 0. Somewhat 
milder but still signihcant restrictions are also placed by on the alternative explanation 
in terms of production of a charm squark off d. Fig. 3 also shows curves corresponding 
to A /21 = 0.04 for tan/3 = 40. The allowed region expands for lower tan/3. Here again, 
other constraints considered in p[ would further limit the available space. In contrast to 
these two case, the process e~^d —>• is not signihcantly constrained by rriu^ as bounds on 
the relevant coupling A /33 are weak due to appearance of Gabibbo angles. It is clear that 
while the neutrino mass bounds may not rule out the SUSY interpretation of HERA events 
they provide by far the strongest constraints which should be systematically included in any 
analysis trying to understand HERA events as the production of a squark. 

All our considerations are based on specihc set of soft terms which are included in MSSM. 
Many of the discussions and in particular the generation of sneutrino VEV in the presence 
of a single trilinear coupling would be true in a more general set up, e.g. in case when soft 
terms arise from the gauge mediated supersymmetry breaking. The main point following 
from our analysis is that the electron neutrino mass provides much stronger constraints on 


8 











R violation and associated phenomenology than has been hitherto realized. 
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0.06 


Fig 1 a ( gi > 0 ) 



Mg (GeV) 

Figure la. FCNC constraints [17] on A '^21 a) m = 200 GeV and h) m = 50 GeV for 
tanjS = 40. Neutrino mass constraints on A '^21 c) m = 50 GeV, tanjS = 15; d) m = 200 
GeV jtanP = 40 and e) m = 50GeV, tanfd = 40. f) Agj^ from neutrino less f3f3 decay g) Agj^ 
from neutrino mass constraints for m = 50 GeV and tanjd = 40. 
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Fig 1 b ( /X < 0 ) 



M2(GeV) 

Figure lb. Same as of Figure la) but with /i < 0 
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Fig 2 a (/i, > 0) 



Figure 2a. Neutrino mass constraints on A '^32 for a) tanf3 = 5 and b) tanjd = 25; on 
for tanP = 5 c) considering only loop contributions and d) loop as well as sneutrino VEV 
contributions are shown for m = 50 GeV. 
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Fig 2 b (// < 0) 



Figure 2b. Same as Figure 2a) but with /i < 0. 
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Figure 3. The horizontal lines correspond to contours of = 85GeV for tan(3 = 40 and 
fi > 0(small dashed), tan(5 = 10 and p, < 0 (dash-dot), tan(3 = 5 and p < 0(dashed). 
Contours of A'^32 = 0.3 for tanf3 = 5 and p < t ) (continuous line), tanf3 = 10 and p < t ) 
(dash-double dot) are also shown. The dotted line corresponds to A 321 = 0.04 for tan(3 = 40 
and p > 0. The region to the left (below) of the contour (the horizontal line) is excluded. 
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